NJL model of homogeneous neutral quark matter: 
Pseudoscalar diquark condensates revisited 
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We use a Nambu-Jona Lasinio type model to investigate the phase diagram of dense quark mat- 
ter under neutron star conditions in mean field approximation. The model contains selfconsistently 
determined quark masses and allows for diquark condensation in the scalar as well as in the pseu- 
doscalar channel. The latter gives rise to the possibility of K° condensation in the CFL phase. 
In agreement with earlier studies we find that this CFLK phase covers large regions of the phase 
diagram and that the predominant part of this phase is fully gapped. We show, however, that 
there exists a region at very low temperatures where the CFLK solutions become gapless, possibly 
indicating an instability towards anisotropic or inhomogeneous phases. The physical significance of 
solutions with pseudoscalar diquark condensates in the 2SC phase is discussed as well. 



I. INTRODUCTION 

Matter at high density and low temperature is ex- 
pected to be a color superconductor. (For reviews see, for 
example, Refs. [B-0])- At very high density, in the limit 
where QCD is weakly coupled, it is possible to investi- 
gate color superconductivity from first principles. From 
these studies [1, Q we know that the ground state of 
three-flavor matter at asymptotically high density is the 
Color-Flavor Locked (CFL) phase [lOj, where quarks of 
all colors and flavors are paired in a homogeneous scalar 
(J p = + ) diquark condensate. Unfortunately, these 
analyses cannot directly be applied to describe the core 
of a neutron star, which is the most likely place for color 
superconducting quark matter to exist in nature. Al- 
though densities of several times nuclear matter density 
can be reached there, this is far from being asymptotic, 
and the weak-coupling expansion breaks down. 

A second difficulty arises from the fact that, at these 
"moderate" densities, the strange quark mass cannot be 
neglected against the chemical potential. Moreover, in 
compact stars the system must be in beta equilibrium 
and electrically and color neutral. As a consequence, 
the Fermi momenta of the different quark flavors are un- 
equal in the unpaired system, and their BCS pairing is 
only favored if the free-energy cost to adjust the Fermi 
surfaces is overcompensated by the pairing energy. It is 
therefore unclear whether quark matter under compact 
star conditions is color superconducting and, if so, which 
pairing pattern is most favored. Besides the CFL phase, 
the most prominent candidate is the 2SC phase (ill Il2|. 
where only up and down quarks are paired. 

The stress imposed by the different Fermi momenta 
of the quarks in a diquark pair reduces the gap in the 
excitation spectrum. Eventually, the gap vanishes com- 
pletely [ljl, and one obtains a "gapless color supercon- 
ducting state" , e.g., the gapless 2SC (g2SC)phase 0OJ] 
or the gapless CFL (gCFL) phase [H, Il7j . However, 
at least at low temperatures, these gapless phases suf- 
fer from chromomagnetic instabilities fla - [20j , indicating 
that they cannot be the true ground state of the sys- 



tem. Possible alternatives are inhomogeneous phases, 
like mixed phases I2l1- |23j or crystalline color supercon- 
ducting phases |24l - [32| . 

Another way the CFL phase can react on the stress 
induced by the strange quark mass is to develop a K° 
condensate. The resulting phase is called CFLK phase. 
Since the K° has strangeness +1 (opposite to a strange 
quark, which has strangeness —1), its condensation effec- 
tively reduces the strange quark content of the system. 
The K° condensation was predicted first within an effec- 
tive field theory approach based on the symmetry break- 
ing pattern of the CFL phase [33- -35] : The CFL conden- 
sates break the original SU (3) color x SU (3) L x SU (3) R 
symmetry of QCD down to a residual SU (3) color+v - sym- 
metry. As a consequence of the spontaneous breaking 
of chiral symmetry, there is an octet of pseudoscalar 
Goldstone bosons, which are massless in the chiral limit. 
Turning on finite quark masses, the Goldstone bosons 
remain rather light [361 ] whereas the strange quark mass 
acts as an effective strangeness chemical potential, which 
eventually leads to Bose condensation of the kaonic 
modes. 

Since the kaon condensate reduces the stress in the di- 
quark pairs, it is expected to shift the occurrence of gap- 
less modes to lower chemical potentials. Within effective 
theories it was predicted that the first mode which be- 
comes gapless is electrically charged [13, HI] . In Ref . [39| 
it was shown that all Meissner masses stay real in this 
phase, but nevertheless a weak instability towards the 
spontaneous generation of Goldstone-boson currents (p- 
wave kaon condensates) was found Ji^, SJ- At lower 
chemical potential a neutral mode could become gap- 
less as well, leading to imaginary Meissner masses and 
stronger instabilities (39j . 

In order to decide, which of these various phases are 
realized under compact star conditions, detailed calcula- 
tions are necessary. Unfortunately, as mentioned in the 
beginning, QCD calculations based on a weak-coupling 
expansion are not applicable at moderate densities. In 
this situation the most promising development is prob- 
ably the Dyson-Schwinger approach, which has been 
successfully extended to describe color superconducting 
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quark matter [42H45| . However, because of the complex- 
ity of this method it is difficult to apply it to the calcu- 
lation of equations of state or phase diagrams. So far, 
the most complete investigations of quark matter under 
compact star conditions have therefore been performed 
within models. Although model calculations have only 
limited predictive power, they can deal as explorative 
studies and may give valuable hints about possible in- 
teresting scenarios or problems. Among these models 
for color superconducting quark matter, the most promi- 
nent ones are models of the Nambu-Jona-Lasinio (NJL) 
type, because they allow for a simultaneous and selfcon- 
sistcnt treatment of the various condensates. Besides the 
diquark condensates this also includes the chiral (quark- 
antiquark) condensates, which determine the dynamical 
quark masses. 

Phase diagrams for homogeneous electrically and color 
neutral quark matter in NJL-type models have been cal- 
culated in Refs. [46|-|49(. Both, diquark and quark an- 
tiquark condensates, were included, but the authors ne- 
glected the possibility of Goldstone boson condensation. 
A rich phase structure was found, including large areas 
of gapless color superconducting phases. Unfortunately, 
the details strongly depend on the model parameters. 

For the description of CFL+Goldstone boson con- 
densed phases in NJL-type models one has to take into 
account pseudoscalar diquark condensates [50|, [HJ- In 
the chiral limit the Goldstone bosons connect the contin- 
uous set of degenerate ground states which are related to 
each other by axial SU (3) transformations, 
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(1) 



where q is a quark field and r a , a = 1 , . . . , 8, are the Gell- 
Mann matrices in flavor space. Under these transforma- 
tions the scalar diquark condensates of the CFL phase 
are partially rotated into pseudoscalar ones. In partic- 
ular, the K° condensation is marked by the formation 
of pseudoscalar condensates which correspond to Eq. ([T]) 
with a = 6, 7, whereas charged pion or charged kaon con- 
densates correspond to a = 1, 2 or a = 4, 5, respectively. 

In Ref. 1531. Warringa extended the NJL-model stud- 
ies of Ref. [46] in this direction. Using the same model 
and the same parameters but additionally allowing for 
pseudoscalar condensates, he found that large parts of 
the CFL and gCFL phases are replaced by the CFLK 
phase. Moreover, Warringa reported that the CFLK 
solutions contain no gapless modes, suggesting that no 
instability problems occur. On the other hand he found 
that, below a certain threshold in the chemical potential, 
the electric charge chemical potential at T = is nonzero 
in this phase. This would mean that there is a nonvan- 
ishing electron density, which is not allowed in the fully 
gapped neutral CFLK phase at zero temperature. 

The main motivation of our present study was to re- 
solve this puz zle. To that end we use the same model 
as in Ref. [52j and carefully analyze the phase diagram 
and the excitation spectra. As we will see, to a large ex- 
tent our results are in good agreement with Ref. 52| . The 



main difference is that we find a region where the CFLK 
phase is gapless. However, this region is restricted to ex- 
tremely low temperatures, T -C fO keV. We will also 
analyze the so-called p2SC phase, a 2SC-like phase with 
pseudoscalar condensates, which was found in Ref. [52| 
but has no obvious physical interpretation. 

The remaining part of this paper is organized as fol- 
lows. In Sec. [IT] we introduce the model and briefly sum- 
marize how the thermodynamic potential is obtained. 
After that, in Sec. IIIIl we discuss the resulting phase 
structure. We begin with a short overview in Sec. IIII Al 
and then perform a detailed analysis of the most inter- 
esting features. Finally, in Sec. IIV1 our results are sum- 
marized. 



II. MODEL 

A. Lagrangian 

We use the same NJL-type model as in Refs. [HI, [k 
The Lagrangian reads 



C = q(ift - m + 7oA)<7 + £gg + C qq + C 



(2) 



where q is a quark field with three color (r,g 7 b) and three 
flavor (u,d,s) degrees of freedom, fa is the diagonal mass 
matrix fa = diag^ (m u , m<2, m s ). 

The diagonal matrix /t contains the quark chemical 
potentials. As motivated in the Introduction, we are 
interested in the electrically and color neutral ground 
state in beta equilibrium. We then have four different 
conserved charges: the total quark number, the electric 
charge and two color charges. According to the cor- 
responding charge densities n, Uq, — n r — n g and 
ng = (n r + n g — 2nb) / \/3, we thus have four independent 
chemical potentials /z, ^3, fig, and /xq. So fi, takes the 
form 



(3) 



with Q = diag / (2/3, -1/3, -1/3) and the Gell-Mann 
matrices A3 and As in color space. 

The other terms in Eq. @ describe the interaction. 
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is a four-point interaction in the quark-antiquark chan- 
nel, while 



,,,< - H [ (qi~f 5 TA^BCq T ) (q T Ci-i5T A \ B q) 

A,B=2,5,7 

T\ f„T, 



(qT A \ B Cq T ) (q T CT A \ B q) ] 



(5) 



is a four-point interaction in the quark-quark channel. 
G and H are dimensionful coupling constants. In our 
notation A a and r a are the Gell-Mann matrices acting 
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The most symmetric phase is the CFL phase, where 



TABLE I: The phases investigated in this article, x marks 
the presence of the condensates while — means that the con- 
densate vanishes. 



in colo r sp ace and in flavor space, respectively, where 
To = ^/2/3 1/. C denotes the charge conjugation matrix 
C = i 7 2 7 - 

Both, Cq q and C qq , consist of a scalar part (first term 
on the r.h.s.) and a pseudoscalar part (second term on 
the r.h.s.) and are invariant under flavor U(3)l x U(3)r 
transformations. In order to break the U(l) axial symme- 
try, we include a six-point ('t Hooft) interaction term [53[ 

C 6 = -K {det/ [q (1 + 75) q] + det/ [q (1 - 75) q]} , (6) 
with coupling constant K. 

B. Phases and condensates 



We introduce the following scalar and pseudoscalar di- 
quark condensates 1 



A<& - -2H(q T C l5 T A \ B q) 
A% = -2H(q T CT A X B q) . 



(7) 
(8) 



Depending on which of these diquark condensates take 
nonzero values, we can distinguish various color super- 
conducting phases. The phases which have been inves- 
tigated in this article are summarized in Table HI For 
instance, in the simplest possible color superconduct- 
ing phase, the 2SC phase, only the A 22 ^ condensate is 
present. In this phase only red and green up and down 
quarks form Cooper pairs while the strange and all blue 
quarks remain unpaired (sec Table ILT]) . In the uSC phase 

[54j A 22 and A 55^ are nonzero. Here red and blue up and 
strange quarks are paired, in addition to the pairs which 
are already present in the 2SC phase. 



1 For the comparison with Ref. 52] it is important to note that 
we use the notation of Re f. 1501 for the diquark condensates, 

^flavor, colon 

olor, flavor • 



while in Ref. 



they are denoted in the opposite 



A 22 \ Agg and A^/ are non-zero, so that all quark fla- 
vors and colors are paired in a scalar condensate. As in 
QCD, the presence of these condensates breaks the orig- 
inal SU (3) color x SU (3) L x SU (3) R of the model down to 
SU(3) 

coior+v Among others, this gives ris e to an octet of 
pseudoscalar Goldstone bosons (see Refs. {55U58I ] for the 
NJL-model description of these Goldstone bosons). The 
condensation of some of these Goldstone bosons is then 
marked by the appearance of nonvanishing pseudoscalar 
condensates. For instance, applying the transformation 
of Eq. (JTJ in the K° channel (a = 6, 7) on the scalar con- 
densates of the CFL phase gives rise to the pseudoscalar 



(p) 

condensates A 25 and A 



$ and Af 2 ', and in the CFLtt 



Similarly, in the CFLK* 

(v) 



phase the condensates A 

phase the condensates A57 and A75 appear. 

In the 2SC phase, there are no pseudoscalar Goldstone 
bosons. Nevertheless, it was found in Ref. 52j that in 
some region of the phase diagram a so-called p2SC phase 
is favored, where a pseudoscalar condensate is present in 
addition to A 22 \ Therefore we will consider this possi- 
bility as well. 

On the other hand, we will not consider other combi- 
nations of scalar diquark condensates, i.e., the dSC phase 



(A 22 } ^ and A$ ^ 0), the sSC phase (A 55 ; ^ and 

A$ ^ 0), the 2SC US phase (A^ ^ 0), or the the 2SC ds 

phase (A77 7^ 0). These possibilities were taken into ac- 
count in Refs. [46[ and [52j, but were not found to be 
realized in the phase diagram. So we do not expect them 
to play a role. 

In addition to the diquark condensates, we also intro- 
duce scalar antiquark-quark condensates 



<t> f 



WQf) 



f e {u,d,s} . 



(9) 



way, 



which are responsible for spontaneous chiral symmetry 
breaking at low temperatures and densities. Since chiral 
symmetry is explicitly broken by the mass matrix m, this 
does not lead to new phases in a strict sense. However, it 
gives rise to dynamical "constituent" quark masses (see 
Eq. p7|) below) which can be considerably larger than 
the bare masses contained in rh. 

Expecting that the condensation of pions and kaons 
outside the color superconducting phase does not oc- 
cur, we neglect condensates of the form (qij^Taq). In 
principle, such condensates should be present in the 
CFL+Goldstone boson phases where they break no addi- 
tional symmetry. However, following Ref. [52}, we neglect 
them there as well, expecting that their effect is small. 
Eventually, this should be checked. 

At the end of this section, we would like to add a few 
more details about the CFLK phase and compare them 
with the more familiar CFL phase. Later on, this will 
be useful for the interpretation of the results. In Fig. Q] 
we show the color-flavor structure of the CFLK pair- 
ing. This pattern follows directly from Tables U and ITT] 
If we switch off the pseudoscalar condensates, we arrive 
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FIG. 1: Pairing pattern in the CFLK phase. 



While all pairs listed in Table [TTI are Q neutral in total, 
we can distinguish two cases: either both quarks carry 
zero Q charge (like in the ru — gd pair) or the quarks 
carry the Q charges ±1 (like in the gu — rd pair). Com- 
paring this with Fig. [1] we see that the two unconnected 
groups in the pairing pattern of the CFLK phase just 
correspond to these two cases. 

C. Thermodynamic potential 



at the CFL pairing pattern. There the quarks can be 
grouped into four subsets, which are not connected to 
each other by the pairing. In the CFLK phase, the 
pseudoscalar condensates give rise to additional connec- 
tions, so that only two unconnected subsets remain. The 
smaller one consists of four pairs, connecting four differ- 
ent quark species, the larger one contains six pairs, which 
connect five quark species. 

In Table |TT] we have also listed the eigenvalues of the 
"rotated" electric charges 



Ah 



(10) 



for both constituents of each pair. For the CFL phase it 
is well known that the ground state is Q neutral |l(t |59| . 
This follows from the fact that all quarks are paired and 
in each pair the sum Q\ + Q2 vanishes. Obviously, the 
same is true for the CFLK phase. An important con- 
sequence is that, at T — 0, in both phases color neu- 
tral quark matter is automatically electrically neutral. 
Hence, no electrons are admitted in neutral matter, i.e., 

=0 (sain. 
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TABLE II: Gap parameters, color-flavor structure of the cor- 
responding diquark pairs and Q charges of the constituents. 
The list is restricted to those condensates which are realized 
in the phase diagrams of Sect. Mil 



In order to decide, which phase is favored at given tem- 
perature and chemical potential, we calculate the ther- 
modynamic potential per volume, 



T 

Q=-]nZ., 



(11) 



where Z is the grand partition function. Besides the 
quarks we consider the presence of noninteracting lep- 
tons, 



(12) 



where fli 



H = Q e + f2 M is the sum of the ideal gas contri- 
butions of electrons and muons. For the quark part fl q 
we restrict ourselves to the mean-field approximation and 
proceed in the same way as described in Refs. [46l.l52l.l62l]. 
Therefore, we only summarize the main steps here and 
refer to these papers for further details. 

Linearizing the interaction around the condensates de- 
fined above and introducing Nambu-Gorkov bispinors, 



1 



1 



V2 V Cf J ' 
we obtain the mean-field Lagrangian 

where 

V =2G (4>l + cf> 2 d + 4> 2 S ) - AKcb u ^ 
+ IH £ [ 

A,B=2,h,7 



(13) 



(14) 



.0) 

AB 



(15) 



is a field-independent potential density, and S 1-1 is the 
inverse dressed propagator, which in momentum space is 
given by 
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It is a 72 x 72 matrix, corresponding to three color, three 
flavor, four Dirac, and two Nambu-Gorkov degrees of 
freedom. M is the diagonal matrix of the constituent 
quark masses with the components 



M a = m a - AGcf)a + 2K(j) b q 



(17) 



where (a, b, c) is any permutation of (u, d, s). 

The quark part of the thermodynamic potential in 
mean-field approximation is then given by 



Q q (T,{ f x i }) = -Tj2 



d 3 p 1 

(2tt) 3 2 



Trln 



S 1 (iuj n ,p) 



T 



+V. 
(18) 

The Tr In is evaluated by determination of the 72 eigen- 
values of Making use of a two-fold spin degeneracy, 
the problem can be reduced to calculating the eigenvalues 
of a 36 x 36 matrix [46| . This matrix can be transformed 
further into block-diagonal form, where the number and 
size of the blocks depend on the pairing pattern. For in- 
stance, in the CFL phase we get one 12 x 12 block and 
six 4x4 blocks, whereas in the CFL + Goldstonc phases 
we have one 20 x 20 block and two 8x8 blocks. The 
different blocks can be attributed to quasiparticle exci- 
tations with different quantum numbers. For instance, 
in the CFLK phase the 20 x 20 block corresponds to 
quasiparticles with Q = whereas the two 8x8 blocks 
correspond to quasiparticles with Q — +1 or Q — — 1, 
respectively. 

The resulting eigenvalues appear in the form iuj n — 
ej(p), where Cj(p) are the energy- independent eigenval- 
ues of the Matrix M. = iuj n — 7°S'~ 1 and correspond to 
the dispersion relations of the quasiparticle modes. This 
allows for an analytic evaluation of the Matsubara sum, 
and we finally get 



18 



d 3 p 



2 T In 1 



V. 



(19) 

Here we have used the fact that the eigenvalues come in 
pairs (ej,—6j), so that we can restrict the sum, e.g., to 
the 18 positive eigenvalues, introducing an extra factor 
of 2. The integral is divergent and will be regularized by 
a sharp three-momentum cutoff A. 

In order to find the most stable solution, we have to 
minimize the thermodynamic potential with respect to 
the various diquark and quark- antiquark condensates, 
i.e., we have to solve the gap equations 



on 

dcf>f 



on 



ODi AB ODi AB 



= 0. 



(20) 



with / £ {u, d, s} and A,B£ {2, 5, 7}. In addition, since 
we are interested in quark matter under compact star 
conditions, our solution must fulfill the neutrality con- 
straints 



= 0, n 3 = - 



dn on 

- — = 0, n s = --r — = 0. 

(21) 



Eqs. (j2TJ)) and (|2"Tj) form a set of coupled nonlinear equa- 
tions which determines the values of the condensates and 
the chemical potentials /iq, /13, and fig for any given 
quark chemical potential /1 and temperature T. For the 
evaluation of the derivatives in these equations we fol- 
low Ref . [HI , making use of the diagonalized form of the 
propagator. Finally, we checked for each solution that 
the off-diagonal color densities vanish as well (cf. Ref. [63[ 
and the discussion in Sec. IIII D|) . 



D. Model parameters 

For the numerical calculations, we adopt the parameter 
values of Ref. [§3, 

m u = m,d — 5.5 MeV, 

m s = 140.7 MeV, 
G = 1.835/A 2 , 
K = 12.36/A 5 , 

A = 602.3 MeV, (22) 

which have been fixed by fitting the vacuum values of 
the pseudoscalar meson masses and the pion decay con- 
stant. The same parameter values have also been used 
in Refs. [46[ and [Hj]. For the diquark coupling strength, 
which cannot be fixed in this way, we follow Ref. (46| and 
take two different values, namely H = 3/4 G, which can 
be motivated by the Fierz transformation of a one-gluon 
exchange, or H — G. 

For the leptonic part, we consider massless electrons 
and muons with mass 105.66 MeV. Here we differ slightly 
from Refs. [ill [52| where a physical electron mass m e = 
511 keV was taken. Obviously this is negligible in most 

As we 



< 



cases, but it could play a role at T, |/xq| 
will sec, this situation is just realized in the vicinity of the 
threshold of the gapless CFLK phase. We will comment 
on this at the end of Sec. IIII Bl 



III. RESULTS 

In this section we discuss the phase structure of the 
model. We begin with a brief overview and then analyze 
the most interesting features in more detail. 



A. Phase diagram 



As discussed in Sec. Ill Cl we obtain the phase diagram 
in the \i — T - plane, by minimizing the thermodynamic 
potential at every point of that plane with respect to the 
various condensates, taking into account the constraints 
of beta equilibrium and local electric and color neutral- 
ity. We thereby consider all phases which are listed in 
Tab. U We recall that, while these phases are distin- 
guished by their contents of nonvanishing diquark con- 
densates, the chiral quark-antiquark condensates 4>f are 
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FIG. 2: The phase diagram of charge and color neutral quark matter for "intermediate diquark coupling", H — 3/4 G. Red 
thick solid lines denote first-order phase transitions, green thin solid lines second-order phase transitions, the blue dotted lines 
mark the (dis)appearance of the gap in the excitation spectrum. 



> 

(D 




320 340 360 380 400 420 

H [MeV] 



440 460 



480 



FIG. 3: The same as Fig. [21 but for "strong diquark coupling", H — G. 



500 



selfconsistently determined as well. Strictly speaking, 
since chiral symmetry is explicitly broken in the model, 
this is not related to additional phases. However, follow- 



ing Ref. [46[, we label regions where 4> u and (fid are large 
by 'xSB' and regions where these condensates are small 
by 'NQ'. Both regimes are not necessarily separated by 
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a phase boundary. 

The resulting phase diagrams are displayed in Figs. [2] 
and El In Fig. [2] the diquark coupling constant was cho- 
sen to be H = 3/4 G. To a large extent our results 
agree qualitatively and quantitatively with those of War- 
ringa, who has calculated the phase diagram within the 
same model for the same parameters [52j j : Compared 
with Ref. 146], where only scalar diquark condensates 
have been taken into account, the main effect of con- 
sidering pseudoscalar condensates is the existence of a 
CFLK phase, which replaces most of the (g)CFL phase 
at temperatures below ss 25 MeV. A second effect, 
which was reported in Ref. J52j, is the emergence of a 
p2SC phase, i.e., a two-flavor color superconductor with 
an additional pseudoscalar condensate, replacing a part 
of the 2SC phase in the region around (i = 460 MeV 
and T = 30 MeV. Although this result is confirmed by 
our calculations (see Fig. [TT] below), we have not indi- 
cated the p2SC phase in the phase diagram, because it 
is most likely a model artifact. This will be discussed in 
SecHlDj 

While most of our results agree with Ref. [52j, there 
are basically two exceptions: First, we find that the 
CFLK phase becomes gapless at extremely low tem- 
peratures and /x < 457.2 MeV, whereas it was claimed 
in Ref. [HJ that it is always gapped. We discuss this 
in more details in Sec. MI Bl The second difference con- 
cerns the phase structure between the normal phase (NQ) 
and the CFLK phase at low temperatures. In contrast 
to the results of Ref. [52[ we do not find a uSC phase 
down to zero temperature. Instead we find a situation 
more like in Ref. [46j, where at temperatures lower than 
~ 7 MeV a gCFL phase is realized. This will be discussed 
in Sec. InTCl 

In Fig. [3] we show the phase diagram for the diquark 
coupling constant H — G. This case was not studied 
in Ref. |52j| . However, it was investigated in Ref. (4^ . 
considering scalar diquark condensates only. By com- 
paring that phase diagram with Fig. 02 we see that the 
effects of the pseudoscalar condensates are qualitatively 
similar to the previous case: A major part of the CFL 
phase is replaced by the CFLK phase, which now per- 
sists up to temperatures of about 35 MeV. At very low 
temperatures and \x < 415.6 MeV the CFLK phase is 
gapless. Numerically, we also find again a p2SC phase, 
which replaces the 2SC phase in some region. However, 
as mentioned before, we believe that this is a model arti- 
fact (see Sec. IIII D[) and do not indicate this phase in the 
phase diagram. 
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FIG. 4: The gap in the excitation spectrum and the electric 
charge chemical potential [iq in the (g)CFLK phase as func- 
tions of the chemical potential for H = 3/4 G and various 
temperatures. 



eral, the gaps are nonzero at higher chemical potentials 
and decrease with decreasing chemical potential. At zero 
temperature (red solid line) it eventually goes to zero at 
(i = 457.2 MeV, so that there is indeed a gapless phase 
below that threshold. 

It turns out, however, that this behavior is extremely 
sensitive to the temperature. At T = 1 MeV (yellow 
dotted line), which is often a very good approximation 
to zero temperature, the gap stays of the order of sev- 
eral MeV and there is no indication for a nearby gapless 
phase. At lower temperatures, the gap first follows the 
T = behavior but then levels off and stays finite. Even 
at T = 0.01 MeV (green dashed line), the gapless phase 
is not reached. This seems to suggest, that the gapless 
CFLK phase only exists at zero temperature. 2 

This picture is corroborated by Fig. [5J where the exci- 
tation gaps are displayed as functions of temperature at 
fixed chemical potential. For /i = 450 MeV the numeri- 
cal values (which are explicitly indicated by the points) 
show an almost perfect linear behavior, so that, within 
numerical accuracy, the CFLK phase becomes gapless 
exactly at T = 0. For comparison, we also show the re- 
sults for fi = 460 MeV. In this case the gap levels off and 
stays finite at T = 0. 

In the lower part of Fig. H] the electric charge chemical 
potential \iq is plotted as a function of ji for the same 
temperatures as the excitations gaps shown in the upper 
part. At zero temperature, as explained in Sect. IIIB1 



B. The gapless CFLK phase 

Both, in Fig. [2J and Fig. El there is a gapless CFLK 
phase at very low temperatures below a certain value of 
fi. For the case H = 3/4 G this is illustrated in Fig. |4] 
In the upper part, the gap in the quasiparticle excita- 
tion spectrum is displayed as a function of /a. In gen- 



2 Note that this has nothing to do with the fact that only at T = 
the transition from a gapped to a gapless phase corresponds to a 
real phase transition with a well defined order parameter. Here 
we identify the gapless regimes by inspecting the gaps in the 
excitation spectrum, which can be done at any temperature. As 
one can see in Figs. [2] and [3] other ga ples s phases defined in this 
way do exist at T > 0, see also Refs. IH, S3, El HI ■ 
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FIG. 5: The gap in the excitations spectrum as function of 
the temperature for H = 3/4 G and two different values of /x. 



fiQ is an order parameter, which vanishes in the gapped 
phase and is nonzero in the gapless phase. This is no 
longer true at T > where fj,Q is always nonzero. Nev- 
ertheless, one might expect that at low enough tempera- 
tures fj,Q is very small in the gapped regime, so that larger 
deviations from zero signal a gapless phase. It turns out, 
however, that this is not the case. For T = 1 MeV, 
0.1 MeV, or 0.01 MeV, [iq exhibits a threshold- like be- 
havior, although we do not enter a gapless regime. When 
we compare these "thresholds" with the upper part of 
Fig. |4j we see that they more or less coincide with the 
positions where the excitation gaps start deviating from 
the T = gaps. In fact, as we will discuss, the mecha- 
nism which keeps the excitation gaps finite at T > is 
closely related to the rise of — fiQ. 

To understand this, we take a more detailed look at the 
low-energy quasiparticle excitation spectra and thereby 
compare the gCFLK phase with the more familiar gCFL 
phase. A similar comparison has been performed earlier 
in Ref. [El]], but restricted to zero temperature and to 
the gapped or the threshold to the gapless regime. Here 
we go deeper into the gapless regime and also discuss 
temperature effects. 

In Fig. [5] the quasiparticle spectra at fi = 450 MeV are 
shown for T = (left) and T = 5 MeV (right). The two 
lower panels correspond to the (g)CFLK phase, while in 
the two upper panels the analogous gCFL solutions are 
shown, which are obtained if the pseudoscalar diquark 
condensates are not taken into account. Since Q is a con- 
served quantum number in both cases, the quasiparticle 
excitations carry a definite Q charge. These Q charges, 
which can be 0, +1, or —1, are indicated in Fig.|5]as well. 

Let us begin with the discussion of the gCFL solutions. 
At T — (upper left panel) there are two gapless modes. 
One of them vanishes at two different momenta, sepa- 
rated by a "blocking region" of about 100 MeV. This 
mode has Q — and can be identified with a 6d-particle 
<?s-hole mode. In addition, there is another gapless mode 
(which has Q = +1 and is a 6u-particle rs-hole mode), 



which seems to vanish just at its minimum. On a much 
smaller scale one finds that this is not exactly true, but 
there is a finite blocking region as well. It is, however, 
very small and therefore not visible in the figure. (For the 
analogous mode in the gCFLK phase, this is discussed 
in Fig.®) 

This behavior is well known and has been explained 
in Ref. [Tjj. It was shown that at T — both modes, 
the bd — gs mode and the bu — rs mode, become gapless 
at the same critical chemical potential but behave rather 
differently when u is further decreased: In the bu — rs 
mode, the pair breaking leads to unpaired bu quarks in 
the blocking region, which carry positive Q charge. Color 
neutral quark matter is therefore positively charged and 
must be neutralized by electrons. This requires /xq < 0, 
which, in turn, reduces the stress in the bu — rs pairs 
and keeps the blocking region very small. For the bd — gs 
mode, on the other hand, there is no such mechanism 
because both quarks have Q — 0. As a consequence, the 
blocking region becomes much wider. 

When we increase the temperature to 5 MeV, as shown 
in the upper right panel of Fig. [51 the bd — gs mode re- 
mains gapless, whereas the bu — rs mode receives a con- 
siderable gap. A detailed numerical analysis reveals that 
this mode behaves in a similar way as the lowest CFLK 
mode in Fig. [5j In particular, it becomes already gapped 
if the temperature is increased only slightly above zero. 
This is naturally explained by Fermi smearing effects. In 
fact, it is a well-known phenomenon that modes which 
are gapless at T = can become gapped above a certain 
temperature [IE|. In the present case, where the block- 
ing region at T = is very small, it is plausible that the 
gapless regime is restricted to very small temperatures. 
On the other hand, since the size of the blocking region 
at T = is not exactly equal to zero, the "critical" tem- 
perature should be nonzero as well. It is, however, too 
small to be determined numerically. 

We now compare these results with the spectra in the 
(g)CFLK phase, which are displayed in the lower part 
of Fig. |51 At T = (left), there is again a marginally 
gapless Q = +1 mode with a tiny blocking region around 
its minimum (cf. Fig. [5] below). The essential difference 
to the gCFL phase is, however, that the gCFLK phase 
contains no second gapless mode, in which Q = quarks 
are paired. The size of the gapless regime is therefore 
determined by the lowest Q = +1 mode. The latter can 
be seen as a "close relative" of the bu — rs in the gCFL 
phase. In the (g)CFLK phase a complication arises due 
to the fact that not only two, but four quark species mix 
in the Q = ±1 quasiparticle modes (cf. Fig. [T}. The re- 
sult of the numerical diagonalization is shown in Fig. [7] 
where the composition of the gapless mode is displayed 
as a function of the three-momentum. In agreement with 
Ref. 5l| , we find that the gapless mode is mainly a mix- 
ture of bu and gu particles with rs holes. The blocking 
of the rs quarks in the gapless regime then leads again to 
a surplus of positive Q charge, which forces the blocking 
region to stay very small. As a consequence, this mode 
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FIG. 6: Quasiparticle dispersion relations in the gCFL phase (top) and in the (g)CFLK phase (bottom) at /i = 450 MeV and 
T — (left) or T = 5 MeV (right); H = 3/4 G. The different line types correspond to modes with different Q charges: Q = +1 
(red solid lines), Q = (green dashed lines), Q = — 1 (blue dash-dotted lines). 



and, hence the CFLK phase is gapless only at very low 
temperatures. 

Our results at T = are in qualitative agreement with 
the effective field theory prediction that the first gap less 
mode in the CFLK phase is electrically charged [37H39j | . 
According to these works, there could be a second phase 
transition at lower chemical potential, where a neutral 
mode becomes gapless as well. In fact, we found such 
solutions in some region of the phase diagram, but they 
were never favored. We will come back to this point in 

sec, mm 

Concerning the electrically charged gapless modes, we 
found some quantitative differences to the effective the- 
ory results. In Refs. [13, HH it was predicted that the 
charged mode becomes gapless at M|/(2/x) w 4/3 A, 
where A is the CFL gap in the chiral limit. Compar- 
ing this with the corresponding threshold in the CFL 
phase, Mf/(2/z) « A, and assuming that the strange 
quark mass is the same in both cases, this means that the 
critical chemical potential for entering the gapless regime 
is 25% lower in the kaon condensed phase. In Ref. 51 1 



it was found that the effect is somewhat weaker in an 
NJL model with A = 25 MeV. It turns out that this ten- 



dency continues with increasing coupling strength: For 
H = 3/4 G we find that both transitions take place 
at approximately the same chemical potential, while for 
H = G the critical chemical potential for entering the 
gapless regime is even shifted upwards by almost 10 MeV 
in the kaon condensed phase. Partially, this can be at- 
tributed to the fact that the selfconsistently calculated 
strange quark mass in the CFLK phase is a few percent 
larger than in the CFL phase, so that the rotation into 
the kaon condensed phase reduces the stress not as much 
as for a constant mass. But even if we take the mass ra- 
tio out, there remains a discrepancy. This indicates that 
there are higher-order effects which have been neglected 
in the effective theory approach, but become important 
at stronger couplings. 

We would like to note that the numerical solution of 
the coupled gap equations (f2TJ|) together with the neutral- 
ity conditions (|2ip turned out to be extremely difficult to 
obtain in the gCFLK regime. As obvious from the be- 
havior of fiQ at T — 0, an extrapolation of the solutions 
from the gapped regime to the gapless regime does not 
work. Formally, this is related to the fact that, whereas 
for the gapped solutions the thermodynamic potential at 



10 



1 








0.9 


■ gu 




0.8 


rs 


\ 

\ 






\ 


0.7 






0.6 






0.5 






0.4 






0.3 












0.2 






0.1 









100 200 300 400 500 600 
p [MeV] 




391.2 391.6 392 392.4 392.8 
p [MeV] 



FIG. 7: Composition of the gapless Q — 1 mode in the CFLK 
phase for H = 3/4 G at fx = 450 MeV and T — (cf. Fig. H 
lower left panel). Shown are the squared coefficients of the 
color-flavor components of the eigenfunctions as functions of 
the three-momentum. 



fixed chemical potentials is a minimum with respect to 
the various condensates, for the gapless solutions it is a 
maximum with respect to A 55 , and (j) s . At the un- 
pairing line, these minima and maxima merge, and the 
solutions correspond to inflection points in these param- 
eters. 3 As a consequence, a small variation of the chem- 
ical potential can easily lead to a situation, where the 
solutions of the gap equations disappear completely. Ob- 
viously, that makes an iterative solution of the problem 
in this regime practically impossible. 

Since the gapless solutions are very close to the un- 
pairing line, these difficulties are hard to avoid in realistic 
cases. We have therefore followed the idea of Ref. 17| and 
introduced additional independent "flavors" of electrons. 
Then, in order to neutralize the extra negative charge 
of these electrons, a larger blocking region is needed, so 
that the gapless solutions move away from the unpairing 
line and the solutions are easier to find. The resulting 
dispersion relation of the gapless mode with 250 electron 
flavors at T = and \i — 450 MeV is indicated by the 
green dashed line in Fig. [8] We see that this mode has 
a blocking region of about 0.5 MeV. For comparison we 
also show the solution for a world without electrons (red 
solid line). In this case, there are no negative particles 
which could neutralize the positive charge of unpaired up 
quarks. Therefore the size of the blocking region must 
vanish, i.e., the solution is exactly located on the unpair- 
ing line. In this case it can be found coming from the 
gapped regime. 

On the scale of Fig. [51 the two cases shown in Fig. [5] 
would be indistinguishable and, hence, indistinguishable 
from the real case with one electron flavor. We also found 



This was also found in Ref. |17H for Alt' in the gCFL phase 



FIG. 8: Dispersion relation of the gapless Q — 1 mode in the 
CFLK phase at fj, = 450 MeV and T = for H = 3/4 G. 
Shown are the neutral solutions in a world with no electrons 
(red solid line) and in a world with 250 independent electron 
species (green dashed line). The solution in the real world 
with one electron species would be very close to the solid 
line, but qualitatively look like the dashed line, with a tiny 
blocking region. 



that the solutions for the gap parameters, masses, and 
chemical potentials obtained with 250 or no electron fla- 
vors differ only little. The calculation with 250 electron 
flavors in Fig. [8] mainly serves as an example to prove that 
gapless CFLK solutions with a finite blocking region in- 
deed exist. Otherwise, of course, the solutions without 
electrons should be closer to the realistic case. For most 
practical purposes, we can therefore well approximate the 
latter by a calculation without electrons. In fact, this is 
what we have done in Figs. 6 and 7, and in Fig. 4 for fiQ 
at T = 0. However, we repeat that the solutions for 250 
electron flavors and, hence, for one electron flavor would 
look the same. 

Finally, we want to comment on the role of a finite elec- 
tron mass m e = 511 keV, which we have neglected so far. 
Since T <C m e in the gCFLK phase, this is only justified 
in the regime where 3> m e . As we can see in Fig. 21 
this condition is fulfilled in the major part of the gCFLK 
phase. In particular, the existence of the gCFLK phase 
is not precluded by the electron mass. Of course, m e 
does have a quantitative effect near the threshold to the 
gapped regime and also on the threshold itself. We recall 
that, because of Q neutrality, color neutral quark matter 
is automatically electrically neutral in the gapped CFLK 
phase at T = 0. Therefore no electrons are admitted in 
the gapped regime, because they would spoil the overall 
neutrality. However, for nonzero m e this does no longer 
force fiQ to vanish, but any \/j,q\ < m e is allowed. As 
a consequence, the system can stay in the gapped phase 
until the unpairing line, which in the relevant regime is 
given by the red solid line in the lower part of Fig. SI 
crosses the line /iq = — m e . This shifts the threshold to 
the gapless phase downwards in /i by about 0.4 MeV. 
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FIG. 9: The free energy of different solutions relative to 
the normal phase, 8fl(T, n) — Q(T,fi) — £Inq{T, [i), for 
H = 3/4 G and T = 0.5 MeV as a function of fi. 
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FIG. 10: The fraction n s /n of strange quarks in different color 
superconducting solutions for H — 3/4 G and T — 0.5 MeV 
as a function of \x. In the boxes it is indicated which phase is 
realized in the phase diagram (Fig. [2)| . 



C. The gCFL window between normal and CFLK 
phase 

We have seen in Figs. [2] and [3] that large parts of 
the CFL or gCFL phase at not too high temperatures 
are replaced by the CFLK phase when pseudoscalar 
condensates are taken into account [52j |. As argued in 
the Introduction, the K° condensate reduces the num- 
ber of strange quarks, which is favorable at "moderate" 
densities, where the strange quark mass cannot be ne- 
glected [33], |34[ . Eventually, at very high chemical poten- 
tials, the CFL phase should be the most favored phase, 
but this happens far beyond the range of applicabil- 
ity of our model, which is limited by the cutoff. From 
this perspective it appears somewhat surprising that for 
H = 3/4 G the gCFL phase survives in a small window 
at the lower-/i end of the CFLK phase and low temper- 
atures, see Fig. [3J This result is also at variance with 
Ref. [52] where a uSC phase was found in this regime. 

Our result is illustrated by Fig. [9j where the free en- 
ergy of different solutions relative to the normal phase, 
Sil = fi — Qnq, at T — 0.5 MeV is displayed against \x 
in the interesting regime. While the normal phase and 
the CFLK phase have the lowest free energy at low or 
high chemical potential, respectively, we see that there 
is indeed a small intermediate regime where the gCFL 
phase is favored. The uSC phase, on the other hand, is 
never favored. 

We have mentioned earlier that for low enough chem- 
ical potentials, effective field theories predict the exis- 
tence of CFLK solutions with neutral gapless modes. In 
the following we will call these solutions 'gCFLK * solu- 
tions' for a better distinction from the "ordinary" CFLK 
solutions with charged gapless modes, as discussed in 
Sec-HIH In fact, at T = 0, the gCFLK * solutions con- 
tain both, neutral and charged gapless modes. However, 
as we have seen for the gCFL phase, the neutral gapless 
mode should be far more stable at finite temperature. 



In Ref. [H] it was found that at T = the gCFLK * 
phase is reached from the ordinary gCFLK phase in a 

(s) 

first-order phase transition, where the condensate Ay 7 
drops substantially. This is not very different from our 
CFLK — > gCFL phase transition, which is also first order 
and accompanied by a strong drop of . In fact, the 
authors of Ref. [39] found that near the phase transition 
the free energy of the gCFLK * phase is very close to the 
free energy of the gCFL phase, so that the result that 
the former is favored is not very robust. 

We have therefore searched for gCFLK * solutions in 
this regime. To that end, we have rotated the gCFL 
solution into the K° direction and used the resulting gap 
parameters as starting values for the iterative solution of 
the coupled gap equations and neutrality conditions. In 
this way, we indeed found gCFLK * solutions, but they 
turned out not to be favored. For T = 0.5 MeV, this is 
shown in Fig. [H] as well. 

In order to shed some light on these results, we com- 
pare the fraction of strange quarks in the various solu- 
tions. This is shown in Fig. [TUJ In the fully gapped CFL 
phase at zero temperature n s /n — 1/3, due to the ab- 
sence of electrons |(5(J. At T — 0.5 MeV, this does no 
longer hold exactly. Nevertheless, as one can see in the 
figure, n s /n is still very close to one third (solid line). 
Below [i = 457.1 MeV the CFL solution becomes gapless 
(dashed line). Near this point the strange quark fraction 
drops very steeply with decreasing chemical potential and 
then approaches smoothly the uSC value (dotted). 

The fraction of strange quarks in the CFLK solution 
is indicated by the dash-dotted line. As expected, the 
K° condensation leads to a reduction of the strangeness 
content as compared to the fully gapped CFL solution. 
However, when going to lower chemical potentials, n s /n 
decreases only mildly and is therefore considerably larger 
than the value in the gCFL solution. From this point of 
view, it is less surprising that the gCFL solution even- 
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FIG. 11: Detail of the phase diagram for H = 3/4 G. In 
contrast to Fig. [2] we have also indicated the region where we 
numerically find the p2SC phase to be favored. 



tually gets favored. Of course, the strangeness content 
is only one aspect. The reduction of strangeness in the 
gCFL phase is related to the breaking of pairs and, thus, 
comes together with a reduction of the pairing energy. As 
a result, the CFLK phase persists to rather low chem- 
ical potentials and the actual phase structure strongly 
depends on details of the interaction. 

For completeness we also show the strangeness frac- 
tion in the gCFLK * solution (dash-double dotted line). 
This solution only exists below ijl = 455.3 MeV, where 
it diverges from the gCFL solution by developing pseu- 
doscalar condensates. At lower fi it approaches the 
gapped CFLK solution, and we have strong indications 
that these solutions eventually join. However, the nu- 
merical analysis becomes rather difficult in this regime 
and, since both solutions are disfavored anyway, we did 
not make much effort to work this out more carefully. 



D. The p2SC phase 

As mentioned in Sec. IIII Al we find numerically that 
the 2SC phase gets partially replaced by a p2SC phase 
when pseudoscalar condensates are taken into account, 
see Fig. [TTJ This was first observed in Ref. [52j . It was 
found there that in those regions where the p2SC phase 
is favored, the free energy gain compared with the 2SC 
solution is very small, less than 1 keV/fm 3 . It was also 
noted that the p2SC solution is connected to the corre- 
sponding 2SC solution via an axial color transformation, 

A 7 

q -> exp(i6»y7 5 )(7, (23) 

which partially rotates the condensate into . In 
this case the condensates in the two phases should be 



related as 




Indeed, the numerical solutions were found to be in very 
good agreement with this relation. We confirm all these 
results. 

Nevertheless, the physical significance of this phase is 
not clear. Unlike the CFLK phase, the p2SC phase 
is not obviously related to the condensation of a Gold- 
stone mode. In fact, the axial color rotation, Eq. (|23l) . is 
neither a symmetry of QCD nor of our NJL-model La- 
grangian, Eq. ([2]), not even in the chiral limit. On the 
other hand, for M u = = and vanishing color chem- 
ical potentials, the mean-field thermodynamic potential 
is invariant under Eq. (|23p if the transformation is re- 
stricted to the non-strange sector and only the diquark 
condensates A 22 ^ and A^ are present (see App. [A"|) . As a 
consequence, there is a spurious pseudoscalar Goldstone 
mode in the 2SC phase at mean field, which is related to 
the spontaneous breaking of this accidental symmetry by 

(s) 

A 22 . This suggests that the p2SC phase is an artifact of 
the mean-field approximation. 

Yet, the accidental symmetry does not explain, why 
the partial rotation of A^ into A^ becomes favored in 
a certain regime, when we turn to the physical case with 
nonvanishing up and down quark masses and color chem- 
ical potentials. In this context, it is quite instructive to 
make a comparison with the case of an ordinary (vector) 
SU(3) color rotation, 

q -> exp(i6 a ^)q. (25) 

It has been observed some time ago that, seemingly, 
this symmetry can be exploited to neutralize the 2SC 
phase without introducing color chemical potentials [65l | . 

(s) 

Thereby A 22 is rotated into a different direction in color 
space, so that equally many red, green, and blue quarks 
are paired in the condensate. As a consequence, the color 
charge densities 713 and ng vanish. Moreover, since no 
color chemical potentials are needed, the resulting state 
has a lower free energy than the standard neutral 2SC 
phase with nonvanishing /Ug. However, as discussed in 
Ref. [63j], this result is not correct if the NJL model is 
taken as a model for QCD: A colored state can of course 
not be neutralized by a color rotation. The color charge 
is just rotated into a different direction, and instead of 
Us one finds nonvanishing off-diagonal color charge den- 
sities n a = (q^Xaq) with a = 1,4,6 [63j. In QCD, these 
color charge densities couple to the corresponding gluon 
fields A° al which then neutralize the system for all color 
components. In NJL-type models, where no gluons are 
present, this must be mimicked by introducing the cor- 
responding (off-diagonal) color chemical potentials. Us- 
ing these color chemical potentials to remove all (diago- 
nal and off-diagonal) color charge densities, the resulting 
state is completely equivalent to the standard 2SC phase, 
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neutralized by fig. In particular it is no longer energeti- 
cally favored. 

Coming back to our original problem, the situation is 
quite similar: When we apply the axial color transforma- 
tion, Eq. (l23|) , to the color non- neutral 2SC phase, the 
ns color charge becomes reduced. However, unlike in the 
case discussed above, it does not completely reappear in 
other components n a . Instead, the transformation gives 
rise to a non vanishing axial color density, 



,.(5) = 



(26) 



To be precise, one finds 



1 



"8 



p2SC 



(1 + 3cos#) n s 



2SC ' 



V3, 



n 3 \^ an = — (1 - cos6»)n 8 | 2SC . 



I P 2SC 
,( 5 ) 



V% . n I 



(27) 



is conserved in the trans- 



Since the sum J2 a ( n a 
formation, the color charge remaining in the vector com- 
ponents gets reduced. Hence, since only the vector color 
densities n a are required to vanish, the rotated solution 
can be neutralized more easily. This favors the p2SC 
phase, whereas on the other hand it is disfavored by the 
nonvanishing up and down quark masses, which break the 
spurious axial color symmetry explicitly. The fact that 
the p2SC phase only appears in a certain region of the 
phase diagram could thus be explained by the interplay 
of these opposite effects. 

As shown in Fig. 1121 we indeed find a nonvanishing 

(5) 

rig in the p2SC phase. Since in QCD there are no "axial 
gluons" , which could couple to these axial color charges, 
there is in principle no need to remove them by intro- 
ducing corresponding "axial color chemical potentials" . 
In this sense the axial color transformation would be a 
legitimate way to reduce the free energy of the system. 
On the other hand, as we have seen before, it is only 
an accidental (approximate) symmetry of the mean-field 
thermodynamic potential. It is therefore unlikely that 
this mechanism works beyond mean field, i.e., the p2SC 
phase is most likely a model artifact. 

Finally, we note that our numerical CFLK solutions 
have a small axial color charge density, too. Although, 
unlike in the p2SC phase, the Goldstone modes which 
condense in the CFLK phase have a well-founded phys- 
ical basis, this could be a model artifact as well. As 
already noticed in Ref . , the solutions of the gap equa- 
tions for the scalar and pseudoscalar diquark condensates 
in the CFLK phase are more similar to an axial color 
rotation than to an axial flavor rotation of the CFL solu- 
tions. On the other hand, both transformations lead to 
the same result in the chiral limit, i.e., the difference is a 
higher-order correction in the quark mass. Hence, even 
if the preference of the axial color rotation was again a 
model artifact, this would be only a small effect on top 
of the real physical origin of the CFLK phase. 
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FIG. 12: The axial color density n ( 6 5) for H = 3/4 G at T = 
25 MeV as a function of fi. The values are nonzero in the 
p2SC phase and smoothly go to zero at the second-order phase 
transition to the 2SC phase at fj, = 451.3 MeV. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have reanalyzed the phase structure of 
homogeneous electrically and color neutral quark matter 
within an NJL-type model. Both, quark-antiquark con- 
densates and diquark condensates have been taken into 
account in order to describe dynamical mass generation 
and color superconductivity on the same footing. The 
main focus was on the role of pseudoscalar diquark con- 
densates, which allow for the possibility of kaon conden- 
sation in the CFL phase. 

A similar investigation was already performed some 
time ago in Ref. [HI . We confirm the main result that the 
CFLK phase replaces major parts of the CFL phase in 
the phenomenologically interesting regime. Numerically, 
we also confirm that there are regions where a so-called 
p2SC phase is favored, i.e., a two-flavor superconducting 
phase with pseudoscalar condensates. We have argued, 
however, that this is most likely a model artifact. 

Our main result is the existence of a gapped CFLK 
phase at small enough chemical potential and very low 
temperature, T -C 10 keV. This is in contrast to Ref. [52j . 
where it was claimed that the entire CFLK phase is com- 
pletely gapped. The smallness of the critical temperature 
of the gapless region is due to the fact that the mode 
contains quarks with nonvanishing Q charge. Like in the 
case of the bu — rs mode in the gCFL phase [l^l, this 
keeps the system close to the unpairing line and, even 
at T = 0, the blocking region stays very small. A small 
increase of the temperature is then sufficient to restore 
the gap by Fermi smearing. 

At T — 0, our results are in qualitative agreement with 
effective field theory, which predicts that the first gapless 
mode the CFLK phase is charged [3714391 ] . However, as 
also observed in Ref. [5l|, there are quantitative differ- 
ences. For large coupling constants we even find that the 
CFLK phase gets gapless at larger chemical potentials 
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than the CFL phase. We also identified CFLK solutions 
with neutral gapless modes, which have been predicted 
within effective field theory as well [37l - l39l |. In our anal- 
ysis, however, these solutions were never favored. 

It is well known that the gapless 2SC and CFL phases 
at low temperatures suffer from chromomagnetic insta- 
bilities, signalling that these phases are not the true 
ground state. For the gCFLK phase with charged gap- 
less modes it was shown in Ref. [39| that all Meissner 
masses are real. Nevertheless the authors found a weak 
instability with respect to the spontaneous formation of a 
Goldstone-boson current, corresponding to a p-wave kaon 
condensate [IS EH- This, and the fact that, at least for 
some parameters, we still found regions where the gCFL 
phase is favored lead us to the conclusion that the phases 
indicated in Figs. [2] and |3] do not always represent the 
true ground state of the system. Most promising candi- 
dates for improvements are inhomogeneous or anisotropic 
phases. This includes the above-mentioned p-wave kaon 
condensates as well a crystalline color superconducting 
phases [HH. 
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Appendix A: Axial color symmetry of the mean-field 
thermodynamic potential in the (p)2SC phase 

In this appendix we show that, for M u = Md = and 
vanishing color chemical potentials, the mean-field ther- 
modynamic potential, Eq. (|18[) . is invariant under the 
axial color transformation, Eq. (|23]). if the transforma- 
tion is restricted to the non-strange sector and only the 
diquark condensates A 2 ^ and A 2 ^ are present. 

If A and A 2 ^ are the only non-vanishing diquark 
condensates, the inverse propagator S _1 , Eq. ((TB|), can 
be put in a block-diagonal form with four 6x6 and six 
2x2 blocks. The 6x6 blocks contain the paired up and 
down quarks. Since all colors participate in the pairing, 
these blocks are larger than the corresponding blocks in 
the 2SC phase, which are 4x4 blocks. The 2x2 blocks 
correspond to the (unpaired) strange quarks. Note that 
for M u = Md = and, hence, <f> u = cf>d — 0, the strange 
sector completely decouples from the non-strange sector 
in the (p)2SC phase. Therefore, since we have restricted 
the axial color transformation to the non-strange sector, 
the 2x2 blocks will not be affected by the transformation. 
Therefore we focus only on the 6x6 blocks. 

Without loss of generality, we may choose the three- 
momentum of the quark to be parallel to the z-axis, i.e., 
P = pe*3 . For M u — Md — and vanishing color chemical 
potentials the 6x6 matrices then take the form 



' Po±M« ~P 
p pa ± fa 



(s ± V 

^ tlfc— d r ~ Ug J 



^25 





A (p) 



1 


^25 





Po T Md 


V 





+ ^22 







A (p) 

^25 


-P 


Po T Md 


















4-aW 
=C ^ A 22 


PO ± Mil 


-p 




V 





+ ^22 





-P 


Po ± M« 



(Al) 



r 



and 



(St-Vr-d.) 1 = [( S ^- dr -n g ) 1 with^Od 

>2) 

Here [l u (/id) is the chemical potential for the up (down) 



quarks. Taking the determinant, we obtain 
-l" 



det 



(s ± V 

y J Ub—d r —Ug J 



(Po ±M«) 2 +P 2 



(O, T Md) ±P) ((Po ± Mu) TP) - (A^) 2 - (A 2 P 5 } ) 2 



Wn2 



((Po T Md) TP) ((Po ± M«) ±p)~ (Ay) 



(^f 
(A3)" 



and the analogous result for det[(S d _ u _ d ) 1 ] with \i n 
and [id interchanged. 
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We see that the determinants and, hence, TrliiS 1 1 = 
lndetS 1 " 1 depend only on the sum (A^ 2 } ) 2 + (A^) 2 , 
which is invariant under the axial color rotation, 
cf. Eq. The same is obviously true for V, Eq. ([15]), 

Consequently, the mean-field thermodynamic potential, 
Eq. (fTSj). is invariant, too. 

Note that it was important that we have restricted our- 
selves to the limit M u = Mj = 0. Otherwise f2 would 



not be invariant. In particular, we would then have to 
transform the condensates 4> u and 4>d a s well. For the 
same reason, we have restricted the transformation to the 
non-strange sector. Alternatively, we could have consid- 
ered the case where the strange quark mass vanishes, too. 
However, this would be a rather poor approximation for 
the real situation in the (p)2SC phase. 
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